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ABSTRACT 

It is recently argued that if the Hawking radiation process is unitary, then a black 
hole’s mass cannot be monotonically decreasing. We examine the time dependent parti¬ 
cle count and negative energy flux in the non-trivial conformal vacuum via the moving 
mirror approach. A new, exactly unitary solution is presented which emits a characteris¬ 
tic above-thermal positive energy burst, a thermal plateau, and negative energy flux. It 
is found that the characteristic positive energy flare and thermal plateau is observed in 
the particle outflow. However, the results of time dependent particle production show no 
overt indication of negative energy flux. Therefore, a black hole’s birth cry is detectable 
by asymptotic observers via particle count, whereas its death gasp is not. 


1 The Puzzle of the Black Hole’s Last Gasp 


The hnal fate of a black hole under Hawking evaporation [1] has long been a puzzle. The 
simplest model of Hawking evaporation, assuming thermal emission, is that the process 
(for an asymptotically flat Schwarzschild black hole, in the units G = c= h = kB = ^) is 
governed by the black body equation (up to a greybody factor). 


dM 

dt 



( 1 ) 


where A ~ denotes the horizon area and T 1/M the Hawking temperature of the 
black hole. In particular, if this equation holds for all time, then it means that the black 
hole mass (for an isolated black hole) is a monotonically decreasing function of time. 
However, we know that if the Hawking process is unitary, then at some point about when 
the black hole loses half of its initial Bekenstein-Hawking entropy, information will start 
to “leak out of the black hole”. This feature is captured by the turn-over of the Page 





curve [2-4], This does not mean that the thermal spectrum and Eq.(l) will necessarily be 
modihed at late times, since the thermal spectrum does not necessitate zero-correlation 
between the radiation [5,6]. Nevertheless, taking into account other effects, it has been 
argued that Hawking radiation is only approximately Planck-shaped over an explicitly 
bounded range of frequencies [6]. 

Despite all these subtleties, the usual thought is that Hawking radiation can only 
decrease the mass of the black hole. In particular, Eq.(l) certainly implies monotonic mass 
loss, and assuming this holds even at late times, the black hole has finite life time of the 
order M^. However, it was recently shown by Bianchi and Smerlak [7,8] that, surprisingly 

— under the reasonable assumption that unitarity is preserved — the mass loss of a black 
hole is not even monotonically decreasing. That is, the mass of an evaporating black 
hole can temporarily increase. The physical picture is as follows: in the usual scenario of 
Hawking evaporation, a positive flux of energy is received at infinity, corresponding to the 
mass loss of the black hole. However, a negative flux of energy is also eventually received 
at infinity^. The puzzle is the following: what exactly causes the black hole to “swell up” 

— the black hole’s last gasp [7]? We have no good physical intuition for such a process, 
however in this work we show that the particle emission process gives no indication that 
the black hole is taking its last gasp. 

In other words, the negative energy flux (NEF) is not characterized by any hint in 
the particle emission. Therefore, an asymptotic observer who is equipped with a particle 
detector and only has access to the particle emission at null infinity, does not detect any 
unusual particle emission, which might otherwise have indicated that the black hole is 
undergoing some unusual process. On the other hand, we show that the thermal flux, and 
even the above-thermal burst (which occurs in our model before the thermal flux) — the 
black hole’s birth cry, is detectable via particle count. 

Due to the difficulty with modeling Hawking evaporation, we shall work in 1-1-1 di¬ 
mensions in natural units, h = c = 1. In (l-l-l)-dimensional spacetime, Newton’s con¬ 
stant G has units to conform to the action principle, rather than the usual 

(3-1-1 )-dimensional units of (This T is of course not to be confused with 

temperature.) We leave G unassigned as we will work in flat spacetime. No unique 
Planck time exists because there is no combination of x, y, z for the dimensional analysis^ 
[G^h^c^] = (7)^ ~ ^ benefit to this unsealed time, is that we are 

allowed to choose time-sampling to be as small as we like without running up against 
Planck scale values. We also employ the moving mirror method to model the evaporating 
black hole spacetime. A moving mirror refers to a free field theory in flat space with the 
boundary condition that the field vanishes on the worldline of a reflecting point (the “mir¬ 
ror”) [13]. An accelerating mirror produces particles via a process called the dynamical 
Casimir effect. Recently, the first experimental evidence for observation of the dynamical 
Casimir effect was reported [14]; see also [15]. 

We are interested in two questions: “Does particle production contain features of the 
energy flux?”; and if so, “Does the particle production contain hints of negative energy 
flux?” 

^See also Figure 12 in [9], in which a moving mirror model was used to investigate the black hole 
firewall issue [10-12]. 

^In H-spacetime dimensions, the Planck time is Tpi = (. Note that the D = 2 case is special. 
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2 An Exact Unitary Solution with Notable Attributes 

We first set up the background in Section (2.1) to explore a unitary solution that has 
features that can be probed both in terms of the particle production and the energy flux. 
We then present the solution and check its consistency in Section (2.2). There are three 
notable attributes: a positive energy burst, a thermal plateau, and a negative energy flux. 


2.1 Some Background: Introducing the Moving Mirror 

The moving mirror model is the most simple example of the dynamical Casimir effect 
[16] [17]. It offers a way to understand black hole evaporation. Consider a moving mirror 
as an external boundary condition in 1+1 dimensions, which does not accelerate forever. 
Consider it starting and ending at time-like past inhnity i~, and time-like future inhnity 

respectively, possessing asymptotically zero acceleration in both the far past and far 
future. As such, an asymptotically inertial mirror will contain no horizons or pathological 
acceleration singularities, as seen by the solid black line in Figures (1). The system, as 
we will see, is unitary. 

As a simple theoretical manifestation of black hole evaporation, the moving mirror 
model, has in practice, been very hard to solve for exact trajectories where the Bogolubov 
transformation may be evaluated, even in the (l+l)-dimensional case. Very few soluble 
cases have been found [17,18] and unitary cases are rare [19,20]. 


t 




Figure 1: Left: In this Penrose diagram, the solid black line indicates the asymptotically inertial trajectory, exhibiting no 
horizon. The grey line is the eternally thermal moving mirror. Right: The asymptotically inertial trajectory with a final 
coasting speed of half the speed of light displayed in the usual spacetime diagram. The dashed lines represent the light 
cone, and the dotted-dashed line shows the asymptote of the mirror trajectory. The trajectory example here is the same 
as in the conformal diagram. For comparison, the grey line indicates the enternally thermal moving mirror [18] [20], which 
contains a horizon coinciding with the light cone. 


The field is the massless scalar of the Klein-Gordon equation O'lfj = 0, where the 
moving mirror is imposed at the z{t) position such that the held is zero, ip\z = 0. The 
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field has modes that can be used to expand it: 


'ip= doj' 


dui'Vui' 


+ 


du [b^r]^ + blr]*^] . 


( 2 ) 


The modes, being orthonormal and complete, take the null coordinate form. 


0^/ = (4W) 




(3) 


= (47ra;)-^/2[g-ia./(.) _ ^4^ 

where the functions p{u) and f{v) are ray-tracing functions related to the trajectory of 
the mirror z{t), see [20]. Expanded modes in terms of the modes themselves. 


/ du T /3uj'uj'Pli^\ ) 

Jo 

poo 

Jo 

is possible from the use of Bogolubov coefficients, 

y Vllj') ■> fJui'ut ) h[j')) 

dehned by the flat-space scalar product, expressed in null coodinates, 

roc ^ POO ^ 

i(j)UyVuj) ■=i du 0 *, duVuj + i dv 0 *, r]^. 


(5) 

( 6 ) 

(7) 

( 8 ) 


The Bogolubov coefficients and [Jui'lo also give the operators a^/ and a[[^, in terms of 
the operators b^^ and bl^: 


a^/ = du [a*,^6^ - f3*,^bl] , 


(9) 


b^= du' 




( 10 ) 


The orthonormality of the modes are expressed via. 


(0a;, 0a;') = “(0)), 0))O =S{u-u'), (0^, 0*,) = 0, 


( 11 ) 


iVu^y Vu) = -(C, Vu) = 6{u-u'), {t]^, T]l,) = 0. (12) 

In the unitary case, the von Neumann entanglement entropy is found from Bianchi- 
Smerlak’s formula^ [8] 

S{u) =-^lnp'{u), (13) 

^Here we set the constant c (the central charge for a conformal field theory) in [8] to unity without 
loss of generality. If the Page curve does turn over, it seems that it would be problematic to fix c [21] if 
the evaporation is not adiabatic. However in our model, the entanglement entropy never decreases. See 
below. 
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and can also be expressed in terms of the mirror trajectory (by using p'{u) = 

S'(t) =-^tanh“^[i:(t)], (14) 

where z{t) is the trajectory motion of the moving mirror, and the dot represents the time 
derivative. In the case we consider below, the faster the mirror moves, the greater the 
entropy. Unitary in this sense means the entropy must achieve a constant value in the 
far past and far future (not necessarily zero). This underscores that eternally thermal 
radiation is not the end state. That is, evolution from a possible initially pure state, to a 
hnal mixed state does not occur. In this context, radiation of energy flux stops and one 
can acheive the initial pure state, preserving unitarity. We will return to this issue in the 
next section. 

The energy flux may also be expressed at the mirror in terms of the entropy as a 
function of time: 

2TiF{t) = cosh^[6S'(t)](6^(t)^ tanh[6S'(t)] + *S(t)), (15) 

which is consistent with the entropy-energy relationship of Bianchi-Smerlak [8] 

27rF{u) = 6S'{u)‘^ + S"{u). (16) 


2.2 The Mirror, the Flux, the Energy and the Beta Coefficients 

To determine whether time dependent particle production contains signatures of energy 
flux, we seek a moving mirror trajectory that has prominent or interesting emission char¬ 
acter. We have found such a dynamic solution. A characteristic feature of this mirror, 
besides the NEF, is a central pulse of positive energy flux. The mirror emits this before 
hnally emitting its NEF and dying back down to zero energy emission. We refer to such 
a pulse as the “birth cry” of the black hole. 

Of course, the features of the energy flux are completely trajectory-dependent; we do 
not claim that this model actually corresponds to any realistic evaporating black hole 
spacetime; for this work we only seek a simple model in which we can study whether 
energy flux (in particular, the NEF) is detectable via particle count. It is nevertheless 
worth pointing out that our dynamic solution is new — we know of no unitary solution 
with a positive energy pulse over the thermal line, like the one exhibited here. 

Starting off asymptotically static, this mirror has no asymptote, ultimately producing 
a hnite total energy^. The motion of the mirror is given by the trajectory 


z(<) = -^ln(e^"«‘ + l), (17) 

where 0 < .^ < 1 is the hnal speed of the mirror as t —)■ cx), and k > 0 is an acceleration 
parameter (not the acceleration, a(t) ^ F). The proper acceleration, a = is found 
to be, of course, non-uniform: 


a{t) = 




\—2^‘^ + (^seclF{K^t) — 2tanh(K,^f)) -|- 4] 


3/2 ■ 


(18) 


^This model therefore respects the observation that evaporating black holes should create finite energy 
flux [22-25]. 
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The negative sign on Eq.(18) indicates the motion is to the left. The acceleration has 
asymptotic behavior such that 

lim a{t) = 0, (19) 

t^±oo 

making this trajectory asymptotically inertial. However, although asymptotically inertial, 
and initially asymptotically static. 


lim z{t) = 0, (20) 

t^ — OO 


this mirror does not approach a future asymptotically static resting state because its 
future asymptotic velocity is 

limi;(t)=^, (21) 

making this trajectory asymptotically coasting. The drifting feature of this mirror means 
it is an exact model for a remnant [26], as anticipated by Wilczek in [27]. For a recent 
review of black hole remnants, see [28]. The entropy of this trajectory using Eq.(14) is 
plotted in Figure (2). 

In this particular remnant scenario, note that the entanglement entropy is monotoni- 
cally increasing, i.e., there is no turn-over of the Page curve (this is not a necessary feature 
of a remnant; see [28]). This means that the radiation is never actually purihed. In the 
black hole evolution context, unitarity is preserved in the sense that the pure state re¬ 
mains pure if we consider both the exterior and the interior of the black hole remnant [26]. 
One may object that such a remnant scenario is problematic due to the well-known in- 
hnite production problem (among other objections), but this is outside the scope of the 
current work. Interested readers are referred to [28] for detailed discussions. Note also 
that our remnant model is not in conflict with the argument for non-monotonic mass loss 
of Bianchi-Smerlak [7], since the latter does not require that the entanglement entropy 
vanishes at the end stage of the evolution; as long as the entropy goes to a constant value, 
Bianchi-Smerlak’s argument follows through. 

The energy flux, in terms of any mirror trajectory, z{t), was hrst derived by Davies 
and Fulling [16]. It is 


F{t) := (T„„) 


'z{z‘^ — 1 ) — 32 ^^ 
127r(i; — l)^(i + 1)^ ’ 


( 22 ) 


where the dots are time derivatives. The energy flux, using the trajectory of Eq.(17), is 
therefore: 


m 


37r [(^ — — 1]^ [(^ -|- -|- 1]^ 


(23) 


and contains three major characteristics: a hyper-thermal positive pulse, a thermal 
plateau and NEF, see Figure (2). In order to investigate particle production, it proves 
useful and possible to write the beta Bogolubov coefficient in terms of the trajectory, z{t) 
(For more details, see Good-Anderson-Evans [20]): 


1 7°° 

(5^^, = -^= / dt [e-*-+*+*--^W(a;+i(t) -a;_)] , 

AttVUJUJ' J-oo 


(24) 
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where a;± := w ± u’. This integral can be evaluated with the trajectory of Eq.(17). The 
solution is in terms of an Euler integral of the hrst kind, i.e. the beta coefficient is a beta 
function; 




a/ 


UJUJ' 


2'kk^{oj — oj') 


B 


oj') i((l + + (1 — 


2< 


2< 


(25) 


The total energy emitted throughout the motion, emitted to an observer at is 


K (^ +l)^tanh ^(0 - ^(2^ + 1) 
dSvr ^ + 1 


(26) 


which we can conhrm from 






doj 


(Tuu{t))il - z)dt. 


(27) 


This verifies that these results are consistent. This also adds weight to the notion that the 
particles carry the energy, in the general sense that the total summation of the energies 
of each quanta is equal to the integral over the energy flux: 


u{Ni^)du = / F{u)du. 


(28) 


3 The Energy Flux: A Thermal Plateau and Two Flares 

A period of thermal emission occurs at extremely high coasting speeds, giving a thermal 
plateau, which is, in the limit ^ 1, located at 


lim Fit) 

t^OO 


dSvr 


(29) 


By construction, this is the same as the constant flux produced by the (eternally thermal) 
Carlitz-Willey trajectory [18], which radiates a thermal Planckian distribution of particles 
for all times. A positive energy flare peaks, in the ^ 1 limit, at 


^(^)peak _ 


4 

3^’ 


(30) 


at time 


^peak _ ^ ^ 

2k 2 


1 

2k 


(31) 


The energy flux that exceeds Eq. (29) is above the thermal line and correspondingly we may 
also refer to it as a “hyper-thermal” flare. This is the black hole’s birth cry in this model. 
Even in the case of hnal coasting speeds where the thermal plateau is unpronounced in 
the energy flux, a hyper-thermal flare can exist. For instance, when ^ = 0.99, the energy 
flux exceeds the value in Eq.(29) but there is only a nascent formation of a thermal 
plateau. For low speeds of there is still a local maximum in the energy flux where the 
positive energy peaks but since this energy does not exceed Eq.(29), we do not call it 
hyper-thermal — the cry is a fainter one. 
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Figure 2: Left: Successive plots of the energy flux, F{t) := (Tuu), Eq.(23), from smallest peak to the largest peak, with 
varying limiting mirror speeds, ^ = 0.7,0.8,0.9,0.99 and ^ = 1 — 0.1^ where x = 3, 4,5,8. Notice the three prominent 
flux characteristics: the positive energy hyper-thermal flare, the thermal plateau, and the NEF. The positive energy hyper- 
thermal pulse and the thermal plateau occur only for a very fast final coasting speed. The NEF is a more robust feature of 
the flux than either the thermal plateau or the hyper-thermal flare (though the local maximum on the curve is also robust). 
The acceleration parameter is set to k = 1. Right: The von Neumann entanglement entropy as a function of time, with 
^ = 0.999987 and k, = 1. One sees that it approaches a constant value, underscoring the final asymptotic coasting speed of 
the mirror and the eventual end of Hawking radiation. There is no turn over of the Page curve. 


For ^ very close to the speed of light = 0.99999999 = 1 — 0.1^, no formal limit is 
taken), the energy flux has a minimal negative value 

, (32) 

^ ^ 3 487r ’ ^ ^ 

which is at the bottom of the valley of the negative energy flux. 

4 The Particle Production 

Particle production has a different dehnition than energy flux [19,29]. However they are 
closely related by the trajectory of the mirror. The information encoded in the trajectory 
provides the link between the energy flux and the particle production. It is already 
known that a mirror can radiate zero, or even negative, energy flux while still producing 
particles [29]. 

The exact solution we considered has this interesting “hyper-thermal” flare, defined 
as that instantaneous energy flux which exceeds the thermal line at k^/AStt. It peaks 
before the thermal plateau at high hnal speeds making it a salient feature of the energy. 
Does this behavior shows up in the particle creation as well? We hnd indeed it does. The 
physical interpretation is that, not only do the particles emitted carry the total energy 
radiated in such a unitary situation but that they carry signatures of the instantaneous 
energy flux emission. This is a stronger evidence than that of the consistency agreement 
of the bare total energy emission from the summation of quanta energy of Eq.(27) alone. 
See Figure (3), in which we constructed “packetized” beta Bogolubov coefficients 



(33) 
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using the orthonormal and complete wave packets of Hawking [1] and calculated the 
particle count, 


{Njn) = / , 


>0 










27 ri(a;i —cjo )n 

e-^ 


(34) 


Particles that arrive at can be investigated by obtaining these wave packet coefficients 
straightforwardly from the coefficients fdojuj' of Eq.(25). With the Bogolubov coefficients 
for the packets in hand, the average number of particles produced for given values of n 
and j is determined. They arrive at in the range of frequencies je ^ u ^ {j + l)e and 
in the range of times (27rn —7r)/e ^ u ^ (27m + 7r)/e. Further details on how these packets 
are constructed from the modes themselves or for non-asymptotically inertial trajectories 
can be found in [20]. 


e=l, .j=2, K=l. f=o.5, e=l, j=2, K=l. 

Njn 



Figure 3: Left: The particle production with final drifting speed very close to the speed of light, ^ = 0.999999999999999 = 
1 — 0.1^^. The acceleration parameter is set to k; = 1. Notice the hyper-thermal pulse of particles where the production 
exceeds the horizontal constant emission line. The horizontal line is the exact solution to particle production with constant 
flux thermal emission [20]. Right: The particle production with final drifting speed at half the speed of light, ^ = 0.5. The 
acceleration parameter is also set to k = 1. In order to compare with the left graph which has a much faster drifting speed, 
notice that at half the speed of light there is a more extreme drop off of particle production in time, i.e. increasing the final 
coasting speed produces a thermal plateau. 


5 Conclusion: Birth Cry vs. Death Gasp 

In this work we have shown that particle production and energy flux, evolved in time, 
exhibit manifest signatures corresponding to the instantaneous radiation. In particular, 
we describe a new, explicit, moving mirror solution which exhibits the thermal plateau, 
as well as the positive energy burst (the black hole’s “birth cry”) before it. Both of 
these features can be detected by asymptotic observers with particle detectors. This 
moving mirror trajectory also gives rise to negative energy flux (NEF), and by the result 
of Bianchi-Smerlak [7,8], also to the non-monotonic mass loss of the corresponding black 
hole. The question is then, whether one could also detect such flux from particle count 
at inhnity. 
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We find no signature of negative energy flux in the particle emission. We do not 
perceive a black hole’s dying gasp via particle emission, in stark contrast to its birth 
cry. The lack of signatures of NEF in the particle production is surprising in light of the 
evidence that the other energy flux characteristics have explicit signatures in the particle 
production. One might expect a rapid decrease or increase in particle emission due to 
the NEF, i.e. a “second flare”; or one might also anticipate zero particle production at 
an instant in time that corresponds to the zero energy flux emission point. We have not 
found these anticipations to have materialized. It should be obvious that what one does 
not expect is negative particle production because the Bogolubov coefficients are complex 
conjugated. 



(35) 


leaving the computed occupation number always positive. Still, the NEF is received at 
inhnity accompanied by real particle production but the character of the NEF is “hidden”. 
Presumably, the NEF is carried by the emitted particles, however the particle emission 
in time gives no hint of NEF (or, correspondingly, an increase in the mass of the black 
hole). This detector, which bleeps, upon a detection, reveals no suspicious patterns that 
the black hole is emitting NEF or swelling up, since the particle production, after its 
central peak, monotonically decreases. Absence of evidence is not evidence of absence, 
however, and therefore we cannot rule out the possibility that the signatures of NEF may 
still be exhibited in the particle production. Evidence might emerge readily in mirror 
trajectories which have salient features of energy flux and achieve asymptotically static 
motion. However, it does seem likely that one could not detect negative energy flux by 
simply looking at particle production alone; instead one may need to use some other 
probes. 

It is also worth noting that in a recent work of Abdolrahimi-Page [21], the negative 
energy flux in the case of an asymptotically flat Schwarzschild black hole is shown to be 
extremely tiny. They found that the mass increase of the black hole is less than 0.09% 
of the energy of a single quantum of the energy of the Hawking temperature of said 
black hole at that time, and is therefore unlikely to be detectable, considering that such a 
signature would in addition be swamped by the noise of quantum fluctuations. However, 
note that in our model the aim is purely to study the detectability of negative energy flux 
at inhnity via particle production evolution signatures, regardless of whether the negative 
energy hux is tiny or not. 
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